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We present numerical results on the static quark-anti-quark grand canonical potential in full QCD at non- 
vanishing temperature (T) and quark chemical potential (/i). Non-zero /j,-s are reached by means of multi- 
parameter reweighting. The dynamical staggered simulations were carried out for nj ■ = 2 + 1 flavors with physical 
quark masses on 4 x 12 3 lattices. 



1. INTRODUCTION 

Forces among infinitely heavy quarks separated 
by distance r are of great significance both at 
T = and when they are surrounded by an inter- 
acting medium with temperature T. The lattice 
formulation of their free energy was successfully 
carried out already 20 years ago pQ. Numerous 
lattice results are available in quenched |2] and 
also in full QCD [514] at T > 0. 

There are recent calculations to determine the 
free energy in channels with definite color trans- 
formation properties. The gauge invariance and 
the interpretation of these singlet /octet channels 
in view of |5I(:>| raised some new, interesting ques- 
tions. 

All the previous results were obtained at van- 
ishing baryonic chemical potential. Lattice sim- 
ulations at n ^ are difficult, because they are 
hindered by the sign problem of the fermion de- 
terminant. Recently, new techniques have been 
proposed to study lattice QCD at non-vanishing 
\i |7I8I9I10| . In this paper we use the multi- 
parameter reweighting to reach non-zero \i values. 
We determine the grand canonical potential of the 
heavy quark-anti-quark system up to (i = 140 
MeV in quark chemical potential on N t = 4 lat- 
tices for physical quark masses. We use the same 
configurations as have been used in 

2. FORMALISM 

The free energy of a QCD system at n = 
with a static quark and anti-quark separated by 
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distance r can be expressed as follows pQ: 
exp(-F qq -(T,r)/T + C) = (L(r)tf(0)) T (1) 

where L(r) = -|tr LJ^I 1 ^(r, r) is the trace of 
the temporal Wilson line, i.e. the Polyakov loop 
and the U^t, r)-s are the link matrices in the time 
direction. C is needed for renormalization, and 
depends only on the lattice spacing (see below). 
F q q is the additional free energy due to the pres- 
ence of the heavy quark-anti-quark pair. 

We can ask for the grand canonical potential 
of the heavy quark system, if we place it into a 
medium with T temperature and /i chemical po- 
tential. One can modify the above formula in an 
obvious way 

exp(-$ OT -(T, M ,r)/T + C) = (L(r)L* (0)) Tlli , (2) 

where <& q q is the additional grand canonical po- 
tential due to the presence of the heavy quark- 
anti-quark. 

Both F q q and $ g g contain the self-energy of 
the static quarks. At T = 0, fi = one usu- 
ally normalizes the potential at a given distance: 
e.g. by demanding that V qq (ro) = 0. Here we 
used the Sommer scale (tq — 0.49 fm). The r in- 
dependent shift, which is applied to the potential 
to satisfy the renormalization condition, is just 
the self-energy at a given lattice spacing. This 
T = shift has to be removed from the finite 
T, fi potentials as well. Using these renormalized 
grand canonical potentials the renormalization of 
the Polyakov-loop can also be done (e.g. |12j). 

In eq. we need the expectation value of the 
Polyakov-loop correlator for non- vanishing fi val- 



1 



2 



ues. For moderate chemical potentials one can 
use the multi-parameter reweighting: one gener- 
ates configurations at fi = and reweights the 
generated ensemble in the (3 and [i parameters, 
simultaneously. The new parameters are usually 
chosen to improve the overlap between the sim- 
ulated and the target configurations, or saying it 
in another way to reduce the systematic errors of 
reweighting to as small as possible. Starting from 
the transition point at fj, = the overlap is max- 
imized along the transition line. However one is 
not forced to follow the transition line during the 
reweighting. The transition ensemble contains 
enough information about both the confined and 
the deconfined phases so a constant T reweighting 
is also possible. (Though this reweighting cannot 
reach as far in /i as the one along the transition 
line.) The reweighting formula for the Polyakov- 
loop correlator is: 



(L(r)Lt(0)) 



<L(r)Lt(0) w (7»> 



T c ,0 



{w(T,n)) Tc ,o 



(3) 



After diagonalizing the transformed fermion ma- 
trix the w weights can be calculated for arbitrary 
H values using the explicit formula of |13| . 

3. RESULTS AT FINITE fi 

We studied QCD with 2+1 flavors of dynamical 
staggered quarks at physical quark masses, that 
is at m Uid = 0.0096 and m s = 0.25 on 4 x 12 3 
lattices. 150 000 configurations were simulated 
at the critical gauge coupling: (3 C = 5.1893. The 
scale was set using the Sommer prescription. The 
Wilson-loops were measured on 24 x 12 3 lattices 
for relatively high quark masses. The scale for 
the physical quark mass is obtained from a chi- 
ral extrapolation: r^/a = 1.77(2). Note that the 
Sommer-scale depends weakly on the quark mass, 
therefore the extrapolation is quite safe. 

Since (3 changes along the critical line, it is nec- 
essary to deal with the problem of renormaliza- 
tion. In order to follow our renormalization pre- 
scription described in sec. 2 we calculated the 
T = potentials not only at (3 C but at further 
two (i values. From the values of the potentials 
at a reference distance (in our case r ) one gets 
a [3 dependent shift, which is used to renormalize 
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Figure 1. <& q q along the T — T c (/i = 0) line for 
fia = 0,0.05,0.10,0.15. Higher /i corresponds to 
lower curve. 
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Figure 2. The <& q q as the function of fi for dis- 
tances (r/a) 2 — 1,2,3,20 along the transition 
line. Larger r corresponds to higher curve. 



the potentials along the critical line. 



3 



On Fig. ^the potential is plotted as the func- 
tion of the distance for various chemical poten- 
tials. The reweighting was done along a constant 
temperature line. For higher \x values the poten- 
tial flattens out at short distances, whereas for 
smaller /z values it reaches its asymptotic value 
only at larger distances. 

Fig. Elshows the potentials as the function of \i. 
The reweighting is done along the transition line. 
Until the critical point (which is around /ia « 
0.18) only slight changes can be observed. After 
this point the errors get enlarged. Since there is 
only a small change in /3, the renormalization has 
only a small effect on the result. 



4. SUMMARY 

We have determined the heavy quark-anti- 
quark grand canonical potentials for non- 
vanishing fi values on N t = 4 lattices. We have 
presented results at fixed T and various [i values 
as a function of r. We have found no significant 
change in the potential along the transition line. 
Around the critical endpoint we have observed a 
rise in the correlations. The renormalization at 
T=0/T ^0 and fi=0/[i ^0 was also discussed. 
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Figure 3. 5{L(r)tf(0)) = (L(r)Lt(0)}„ - 
(L(r)ZJ (0)) M =o, i.e. the difference in the corre- 
lator to the fi = case as the function of /i. The 
distance is {r/a) 2 = 13. 

It can be instructive to examine the change in 
the Polyakov-loop correlator along the transition 
line. Fig. [3] shows the difference in the correlator 
to the n = case for a given distance. The corre- 
lator starts to grow intensively before the second 
order endpoint. Unfortunately, at the same time 
the statistical errors get enlarged. 



REFERENCES 

1. L. D. McLerran, B. Svetitsky, Phys. Rev. D 
24 (1981) 450-460 

2. O. Kaczmarek et al, Phys. Rev. D 62 (2000) 
034021 

3. O. Kaczmarek et al., Prog. Theor. Phys. 
Suppl. 153 (2004) 287-294 

4. P. Petreczky, K. Petrov, 
[arXiv:hep-lat /0405009| 

5. V. A. Belavin et al., Phys. Lett. B 579 (2004) 
109-112 

6. O. Philipsen, [arXiv:hep-lat /0407042| 

7. Z. Fodor, S. D. Katz, Phys. Lett. B 534 

(2002) 87-92 

8. C. R. Allton et al., Phys. Rev. D 66 (2002) 
074507 

9. P. de Forcrand, O. Philipsen, Nucl. Phys. B 
642 (2002) 290-306 

10. M. D'Elia, M. P. Lombardo, Phys. Rev. D 67 

(2003) 014505 

11. Z. Fodor, S. D. Katz, JHEP 0404 (2004) 050 

12. O. Kaczmarek, F. Karsch, P. Petreczky and 
F. Zantow, Phys. Lett. B 543 (2002) 41 

13. Z. Fodor, S. D. Katz, JHEP 0203 (2002) 014 



